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ABSTRACT 



Surveys based on the Sunyaev-Zel'dovich (SZ) effect provide a fresh view of the galaxy cluster population, one that is complementary 
to X-ray surveys. To better understand the relation between these two kinds of survey, we construct an empirical cluster model using 
scaling relations constrained by current X-ray and SZ data. We apply our model to predict the X-ray properties of the Planck SZ 
Cluster Catalog (PCC) and compare them to existing X-ray cluster catalogs. We find that Planck should significantly extend the depth 
of the previous all-sky cluster survey, performed in the early 1990s by the ROSAT satellite, and should be particularly effective at 
finding hot, massive clusters (T > 6 keV) out to redshift unity. These are rare objects, and our findings suggest that Planck could 
increase the observational sample at z > 0.6 by an order of magnitude. This would open the way for detailed studies of massive 
clusters out to these higher redshifts. Specifically, we find that the majority of newly-detected Planck clusters should have X-ray 
fluxes lO"''' ergs/s/cm^ < fx[0.5 - 2 keV] < 10"'^ ergs/s/cm', i.e., distributed over the decade in flux just below the ROSAT All Sky 
Survey limit. This is sufficiently bright for extensive X-ray follow-up campaigns. Once Planck finds these objects, XMM-Newton and 
Chandra could measure temperatures to 10% for a sample of ~ 100 clusters in the range 0.5 < z < 1, a valuable increase in the number 
of massive clusters studied over this range. 

Key words. Cosmology: cosmic background radiation; Cosmology: observations; Galaxies: clusters: general; Galaxies: clusters: 
intracluster medium; X-rays: galaxies: clusters 
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The Sunya ev-Zerdovich (SZ) effect dSunvaev & ZeldovichI 
1 19701 Il972h offers a promising new way of studying galaxy 
clusters, one that is complementary to the more tradi- 
tional methods based on X-ray and optical/IR observa- 
tions. A suite of dedicated instruments is already produc- 
ing high-quahty SZ measurements of previously known clus- 
ters ("Udomprasert et al.'2004';'Muchovei e t al.l2007l lZwart et al.' 
,2008; Basuetal. 2010; Halverson et al.l 120091: iHincks et al. 
[20101; lPlaggeetal.1120101; IWu et all l2009i) . and the long antic- 
ipated era of SZ cluster surveys has arrived: the South Pole 
Telescope (SPT), Atacama Cosmology Telescope (ACT) and 
the Planck satellite are all discovering clusters solely through 
their SZ signal dStaniszew ski et al. 2009; V anderlinde et al . 
l2010HMenanteau et al.ll2010 ; Planck Collabora tion et al.ll201 lal) 
In particular, the Planck satellite^ is performing an all-sky 
SZ survey and the Planck consortium has recently published 
an early list of clusters from 1 months of observations 
dPlanck Collaboration et al.l 1201 Ibl) . The Planck survey is the 
first all-sky cluster survey since the ROSAT All-Sky Survey 
(RASS) of the early 1990s. The ROSAT dataset produced sev- 
eral reference catalogs either directly based on the RASS, among 
which NORAS (Northern ROSAT All-Sky galaxy cluster survey, 
iBohri nger e t al. 2000) . REF LEX (ROSAT-ESO Flu x Limited X- 
ray Galaxy Cluster Survey, IBohringer etaLll2004l) and MACS 
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(Massive Cluster Survey, Ebeling et al. 2001), or serendipi- 
tous catalogs, such as the 400sd (400 Square Degree survey, 
iBurenin et id]l2007t) . the SHARC surveys ( Serendipitous High- 
Redshift Archival ROSAT Cluster survey, iRomer e t al."2000l: 
Burke et al. 2003) or the WARPS surveys (Wide Angle ROS AT 
Pointed Survey, iPerknan et all 120021; [Horner et all l2008l) . A 
comprehensive descr iption of these surveys can be found in 
Piffaretti et al.l (l2010l) . The forthcoming SZ surveys will find 
many new clusters in addition to numerous objects akeady 
known from X-ray (and optical/IR) studies. 

It is important to understand the relationship between SZ and 
X-ray cluster surveys for a number of reasons. Firstly, compari- 
son of the two kinds of survey will clarify the nature of the se- 
lection functions of both types of catalog. Secondly, comparing 
X-ray and SZ catalogs will yield important new scientific results 
on cluster physics, including information on the SZ-mass rela- 
tion crucial for cosmological applications of cluster evolution. 
Finally, it will help dimension follow-up X-ray observations of 
newly discovered SZ clusters, most notably those at high red- 
shift. 

With these objectives in mind, we develop an empirical 
model for X-ray and SZ cluster signals based on observed in- 
tracluster medium (ICM) scaling relations and use it to pre- 
dict the expected X-ray properties of SZ-detected clusters. In 
the present paper, we focus specifically on the expected Planck 
Cluster Catalog, which we refer to as the EPCC to distinguish 
it from the actual future. We quantify the expected overlap be- 
tween the RASS and the EPCC, and show how the latter sig- 
nificantly extends the RASS cluster catalogs from z < 0.3 to 
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redshift unity. We find that newly-discovered EPCC clusters at 
z > 0.5 are hot. X-ray luminous systems with apparent X-ray 
fluxes falling in the decade just below the RASS sensitivity: 
10-'^ ergs/s/cm^ < fx[Q.5 - 2 keV] < lO^'^ ergs/s/cml 

This result has important consequences for follow-up obser- 
vations, implying that large numbers of Planck clusters could 
be studied in detail with dedicated programs on XMM-Newton 
and Chandra. As an example, we show that with 25-50 ksec ex- 
posures, XMM-Newton could measure the temperature of EPCC 
clusters to 10% out to z = 1 . A 5 Msec program would then yield 
~ 100 massive clusters with measured temperatures and, in the 
best cases, gas mass profiles in the relatively unexplored red- 
shift range 0.5 < z < 1 . This follow-up would help calibrate the 
SZ-mass relation in this redshift range and extend the reach of 
clust er gas ma ss fraction measurements and dark energy studies 
(e.g.. lGalli et al. (2012)). Such a project falls naturally under the 
category of Very Large Programmes envisaged as legacy projects 
fox XMM-Newton after 2010. 

We begin in the next section with a presentation of our em- 
pirical cluster model and its observational basis. We then discuss 
the EPCC, focussing on expected detection limits and the result- 
ing catalog by using both analytical calculations and detailed 
simulations of Planck observations. Given the expected content 
of the Planck catalog, we apply our model to predict its X-ray 
properties and compare to the RASS. Finally, we examine the 
ability of XMM-Newton to follow-up in detail a large number of 
newly-discovered PZfl«cA: clusters, before concluding. 

We emphasize that none of the results presented here was 
derived from Planck data. The characteristics considered here 
for the Planck survey all correspond to the pre-launch expecta- 
tions, as defined in the Planck Blue Book ( The Planck Consortia 
l2005h . For that reason, the EPCC discussed here should be con- 
sidered as representative of a Planck-\\k& cluster catalog. 

2. Cluster Model 

We seek a simple, empirical model relating quantities that are 
directly observable to each other and to cluster mass and redshift. 
This is sufficient to our purpose of establishing a baseline model 
in order to: 

1. Robustly predict the expected X-ray properties of SZ- 
detected clusters by direct extrapolation of the current ob- 
servational situation; 

2. Provide a reference for the interpretation of combined SZ 
and X-ray studies; deviations from the model predictions 
will provide clues to missing physics, providing important 
feedback to more detailed theoretical modeling; 

3. Better understand the validity of crucial scaling laws needed 
to use clusters as a cosmological probe, e.g., the cluster 
counts. 

Although the model is completely general, in the present paper 
we focus its application on the EPCC, leaving other SZ surveys 
to future work. 



Clusters are essentially dark matter halos hosting hot gas, 
the X-ray emitting intracluster medium (ICM), and galaxies. To 
good approximation, most properties are primarily functions of 
cluster mass and redshift, at least when averaged over the entire 
cluster population. That this is the case is perhaps not too surpris- 
ing given that the formation of their dark matter halos is driven 
by gravity alone. The absence of a preferred scale in gravity then 
implies that relations between a halo property, e.g., virial tem- 
perature, and mass and redshift should be power laws. In fact, it 
implies that expo nents of these power-law relations should have 
particular values (lKaiseilll986l) . Deviations from these average 
relations indicate the importance of other factors, such as accre- 
tion history or shape of the initial density peaks (iGao & WTiitd 
120071: iDalal et al]|200^ . Massive systems like clusters are rare 
and isolated, and present a more homogeneous population than 
lower mass systems, such as galaxies. This is precisely one of 
the reasons they are good cosmological probes. 

iNavarro et alj (1 1997b . hereafter NEW, identified a universal 
dark matter density profile for halos from numerical simulations 
depending on two parameters that can be directly related to clus- 
ter mass and redshift. The NEW profile represents the average 
dark matter distribution for halos, while of course individual ob- 
jects scatter about this average. We adopt the NEW mass distri- 
bution for our model clusters, ignoring any scatter. This is ade- 
quate for discussing average properties of the cluster population, 
but should be kept in mind when dispersion of individual objects 
is under consideration. 

The cluster gas is subject to physics other than just gravity, 
such as cooling and heating by feedback from member galaxies. 
We could therefore expect the relations between gas properties 
and mass and redshift to differ from the self-similar power laws 
applicable to the dark matter halos. In the more massive clus- 
ters, however, gravity dominates the overall energetics and most 
of the gas scaling relations do not greatly deviate from their self- 
similar predictions. Deviations tend, rather, to appear in lower 
mass systems {Jx < 3 keV). These deviations are important 
clues to the gas physics and hints to important processes driv- 
ing galaxy formation. 

The key ingredients of our model are: 

1 . A fiducial cosmology, which we take as the WMAP-7 best- 
fit model (Larson et al. 2011); 

2. Tx - M500 relation dArnaud et al.l 120051; IVikhlininetal] 
l2006h ; 

3. Ly - Msnn relation (iPratt et alJl2009h ; 

4. Y-M 5Q0 relation and profile for the SZ signal dArnaud et al.l 
l20Toh ; ^^_^ 

5. A h alo mass func t ion — we use both the lJenkins et al.l (1200 ih 
and [Tinker et al.l (l2008l) mass functions. Eor the former, we 
need a relation between M500 and the halo mass (detailed 
below). As we shall see, the two mass functions yield very 
similar predictions once normalized to the local cluster abun- 
dance. 



2.1. Approach 

We build the model on observed scaling relations between X-ray 
or SZ observables and cluster mass and redshift, which are the 
fundamental cluster descriptors. The properties of the dark mat- 
ter halos hosting clusters are taken from numerical simulations 
and form the scaffolding on which we construct the model. This 
is a standard approach, but we invest some time in its description 
to highlight its strengths and limitations. 



As discussed below in Sec. 12.61 the quantity M^^ refers to the 
cluster mass within the region where the mean mass density is 
Ac times the critical density at the redshift of the cluster, pdz)', 
the radius of this region is consequently referred to as ta . These 
quantities are then related by: Ma^. = 3^'"^ Pciz)l!^c- The param- 
eters of the mass function fit by Tinker et al. (2008) are given 
as a function of the mean overdensities A„, = -^^-j- ; as a conse- 
quence, we translate our critical overdensities A^ to A,„ in order 
to use this mass function in the appropriate terms. 
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We adopt the standard flat ACDM WMAP cosmology with 
Qm = 0.262 = 1 - Qa and Ho = 7L4_km/s/Mpc, taken 
from the WMAP seven-year analysis dLarson et al. 201 1|), and 
define /170 = Hq/{7Q km/s/Mpc). We explore, however, dif- 
ferent values for erg, incl uding the WMAP-7 value {os - 
0.801 dKomatsu et al.ll201ll) ) as well as our best fit to the counts 
in the 400 Square Degree Survey (400sd) ("Bu reninet ani2007h 
with the observed Lx - M50 relati on 02.3b . The latter, com- 
bined with the Jenkin set al.l ( 1200 Ih mass function, leads to a 
slightly lower erg, although still consiste nt to within 1.5cr with 
the WMAP-only value; the [Tinker et all (|6008) mass function, 
on the other hand, leads to fully consistent results. 

We now describe each of these ingredients in turn. 



2.2. Mass-Temperature Relation 

Clusters are usefully described by a global X-ray temperature, 
although they are certainly not isothermal, displaying a large 
temperature scatter in the c ore region and a declining tempera- 
ture profile in the outskirts ( Pratt etal.l2007h . X-ray observations 
of galaxy clusters provide both gas temperature measurements 
and total mass estimates, the latter t hrough the as s umpti on that 
the gas is in hydrostatic equilibrium. ' Arnaud et aL ( 2005b. using 
XMM-Newton observations of 10 clusters, and IVikhlinin et aO 
(1200 6). using Chandra observations of 13 clusters, indepen- 
dently measured the X-ray mass-temperature relation for re- 
laxed systems at z = 0. Following these authors, we parameterize 
the mass-temperature relation as 



E{z)Msoa = Ms ( 



T \'-*' 



5keV 



) 



h-1'0 



(1) 



where E{z) = ^ - [O.ui'^ +1}^ + ^a]'^^ and where we take for 

our fidicual values M5 = 4.1xlO''*Mo anda^r = 1.5, consistent 
with both of these analyses as well as with self-similar scaling 
in cluster mass. 

2.3. Luminosity-Mass Relation 

This is one of the most important relations because it links the 
mass and redshift to the basic X-ray observable, the luminosity. 
We will, in fact, use several Lx - M relations depending on the 
situation, e.g., energy band considered, extent of the region in 
whic h the luminosity is measured, etc. They all, however, come 
from iPratt et al.l (l2009l) and w ere derived from the R EXCESS 
sample dBohringer et al.ll2007h . We write the relation as follows: 



EizT^I^L = C 



M 



2 X W%^^M^ 



hll 



(2) 



with C and aiM given in Table[T] This relation shows a clear de- 
viation from the self-similar model, where the slope olm should 
be 4/3. In our model, this deviation is explained by a temperature 
dependence of cluster gas mass fra c tion consistent with a variety 
of obs ervations (e.g.. iMohr et alj dl999l) : iNeumann & ArnaudI 
d2001h V 

2.4. SZ Signal 

The SZ signal is measured as a surface brightness change to- 
wards a cluster relative to the mean sky brightness (the cosmic 
microwave background, henceforth CMB): diy = yjv where jy 
is a universal function of the observation frequency only and 
the amplitude is given by the Compton-y parameter y(h) - 



J dl (kTe) / im^c^)necrj , an integral of the electron pressure P - 
riekTe along the line-of-sight in the direction given by the unit 
vector «; the e subscript refers to electron q uantities and ctj is the 
Thomson cross section dBirkinshawll 19991) . Because the spectral 
signature is universal, we can define the SZ "flux" in terms of 
the integrated Compton parameter Y - J dQ.y(h), i.e., integrated 
over the cluster image. This quantity has no physical units and 
we express it in arcmin^. 

For our SZ model, we employ the gas press ure distribu- 
tion P(r) recently deduce d from X-ray observations ("Nag ai et al] 
2007; Ar naud et 301201 Oh and specified as a mo dified NFW pro- 
file. We use the results of Arnaud et al.l d2010), who fit a self- 
similar pressure profile to the REXCESS sample to obtain: 



P(r) 

P5m 



Po 



x^il + xP^)<P-y)i" 



(3) 



where x - r/r^ with r, - r^oo/csoo and C500 - 1 . 156, a = 1 .0620, 
/3 = 5.4807, r = 0.3292. The normahzation Pq = S.UOh^^y The 
quantity P500 expresses the scaling and is given its self-similar 
value: 



P500 = 1.65xl0-^£(z)**/^ 



M500 



3xl0i4/i-jiMay 



2/3 



'70 



keV cm"^ (4) 



Assuming 0.3 solar metallicity for the ICM (the result 
changes very little with gas metallicity) and integrating out to 
a projected distance of 5 x rsm, we calculate the Y - M500 rela- 
tion: 



Y [arcmin^] = Y* 



500 



M500 



3 X lO'^o'Mc 



5/3 



Eiz) 



2/3 



Dangjz) '^ 

500Mpcy 



(5) 



2.87 X 10 arcmin and Da„g is the angular 



where Y*^ 

diameter distance. This pressure profile converges so that the 

choice of outer radius does not greatly a ffec t the r esults. 

As shown by iMelin et al] d2011h and in 
iPlanck Collaboration et al.l d2011dld) . this model provides 
an excellent match to direct SZ observations. 



2.5. Redshift Evolution 

The redshift evolution of these scaling laws remains, unfortu- 
nately, quite poorly constrained, mostly as a consequence of the 
lack of high-redshift data. Recent studies are, however, consis- 
tent with self-similar evolution of the M^po - Ty and Ly - Tx 
relations (iKotov & VikhUninll2006t iPacaud et ai1l2007l) . And in 
the absence of any constraints on the evolution of the M500 - Y 
relation, we adopt self-similar scaling in this case as well. We 
then have: 



M500 I 
oc E (z) 

L 



l^aLM 



oc E'^I\z) 



7Vi 



-2 



M. 



<x E''\z)D-J(z) 



(6) 
(7) 
(8) 



500 



as written explicitly in our expressions above. 



2.6. Mass Conversion 

Numerous definitions of cluster mass are used in the literature, 
the reason being that there is no unique definition of cluster ex- 
tent. As a consequence, our model constraints and input relations 
employ different mass definitions. 
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Table 1. Luminosity-Mass relations used in the model, expressed in the form E{zy''^^L - C I ,y_i I /z^q . These relations are 

given in iPratt et alJ (l2009l) . from analysis of the REXCESS dBohringer et al.ll2007h data. All relations are corrected for Malmquist 
bias. The last column gives the logrithmic dispersion of the luminosity relation. 



Application 



Energy band [keV] Region [rsoo] C [10 erg s ] 



TTi 



*-^lnL, intrinsic 



400sd distribution [0.5-2] 0-1.0 0.48 ± 0.04 1.85+0.14 0.414± 0.071 

REFLEX distribution [0.1-2.4] 0-1.0 0.78 ± 0.07 1.85+0.14 0.412± 0.071 

EPCC properties for XMM [0.5-2] 0.15-1.0 0.38 + 0.02 1.53+0.05 0.174± 0.044 



We can basically divide them into three types: masses used 
in theoretical relations (e.g., M,.,,., the mass contained within r,,,>, 
the virial radius), masses used in (most of the) observed rela- 
tions (e.g., Msoo) and masses employed when analyzing numeri- 
cal simulat ions (e. g., masses estimated by the friends-of-friends 
technique (iDavis e t al. 1985), hereafter Mf„f). The first two are 
defined via a chosen value for the density contrast Ac, defined 
such that Ma^ is the mass over a region within which the mean 
mass density is Ac times the critical density at the redshift of the 
cluster, p - A^^pciz) - Ac 3//^(z)/87rG . For the virial mass, 
M^iir, the density contrast Ac varies with redshift, while in the 
second case it is fixed, e.g., at Ac = 500; this is considered a 
good compromise between the inner regions affected by core 
physics and the potentially unrelaxed regions near the virial ra- 
dius (Ac - 180 for an Einstein-de Sitter cosmology). 

In contrast, the friends-of-friends technique establishes the 
perimeter of a cluster using a threshold in the distance between 
closest neighbouring particles. This defines an isodensity con- 
tour around the cluster, without any prior on the dark matter dis- 
tribution in the cluster itself. The isodensity contour refers to 
the background density in the simulation, rather than the critical 
density, at the cluster redshift. 

Eventually, everything needs to be expressed in terms of the 
same ma ss, say Mfpf, since this is the mass used in the mass 
function (I Jenkins et al.ll200lh . To convert between the various 
mass definitions, we require a dark matter profile. To this end, we 
adopt an NEW dark matter profile with concentration C200 - 4.3 
(This concentration value corresponds to the error-weighted av- 
erage for the clusters used by Vikhlinin et al. (2006) to derive 
their M500 - T relation). Note that the conversion factor between 
mass definitions depends on redshift as well as on the under- 
lying cosmology. Simulations furthermore indicate that hydro- 
static mass estimates are biased low due to non-thermalized bulk 
gas motions (Piffaretti & Valdarnini 2 008; Arnaud et al. 2010t). 
In our mass conversion we account for a 15% bias between the 
hydrostatic mass used in the scaling laws (e.g. M500) and M/of. 

2.7. Independent Validation 

As a check of the model, we compare its predictions to other, 
independent observational constraints, starting with a look at the 
X-ray c luster counts. Using the Lx - M relation (see Table [TJ 
and the I Jenkins et al.l (l200 ^ |) mass function, we fit to the 400sd 
survey counts (Burenin et al.ll200 7) to find os = 0.763 + 0.008, 
keeping all other paramters fixed at their fiducial values given 
above. The result is shown in Eigure[T] Using the iTinker et al.l 
(120081) mass function, we obtain as = 0.801 + 0.009. Both 
values agree well (within 1.5cr in the worst case) with the 
WMAP-7 result of os = 0.801 + 0.030, with q-jj ^ 0.773 + 
0.025 , value derived from the SPT power spectrum ( Lueker et al.l 
I2OIOI) . and lie within the range allowed by most recent measure- 
ments dTvtler et al.l2004tlCole et al.l2005l:lReichardt et al.l2009t 
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Fig. 1. Observed r edshift distribution ( continuous black line) in 
the 400sd survey dBurenin et al.ll2007|). compared t o the model 
predictions for cr^ - 0.763 with the iJenkins et al.l ( 200 lb mass 
function (dashed blue line) and erg = 0.801 with the lTinker et al.l 
(12008 ) mass function (dot-dashed red line). 



iJuszkiewicz et al.ll2010b . Eurthermore, IVikhlinin etal] (l2009bl) 
estimated the relation between erg and ^m using clusters from 
the 400sd survey and the RASS. They found: 



0"g 



025 



-0.47 



0.813 + 0.013, 



(9) 



leading to erg = 0.795 + 0.013 when using Q.yi - 0.262 as we 
did for the rest of this study, in full agreement with our estimate. 
This is a non-trivial and important consistency check. 

The inclusion of the 15% bias in mass proves indispensable, 
as our best estimates for erg, when omitting this bias, drop to 
erg ^ .72 2 + 0.008 and 0-8 = .759 + 0.008 for the J enkins et al.l 
(1200 ih and [Tinker e t al.l d2008l) mass functions, respectively. 

To follow the remaining uncertainty associated with erg and 
the mass function in the rest of this paper, all estimates will be 
derived using erg - [0.75, 0.80, 0.85] for both mass functions. 

Secondly, the predicted local X-ray Luminosity Function 
hereafter XL E) is ful ly consistent with the measured XLE 
lEbeling et al. 1997; Bo hringer et al.l2002h . as shown on the left- 
hand side of Eig.|2] This is in particular true at high luminosities, 
i.e., for those clusters of greatest interest to our present study. 

Next, the mod el reproduces the reds hift distribution of the 
REELEX survey (iBohringer et al.l 2004h . a flux-limited survey 
(with fdet = 3 X lO-'^ergs/s/cm^ in the [0.1-2.4]-keV band) 
covering a total area of 4.24 ster in the Southern Hemisphere. 
The survey contains 447 optically confirmed clusters and claims 
to be at least 90% complete. Eor an equivalent survey, our 
model predicts a total of 507 clusters for erg - 0.763 (using 
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Fig. 2. Consistency checks of the cluster model. Left: L ocal predicted XLF f or the two mass functions (along with fitted 0-^ val- 
ues), c ompared to obse rvations by Ebelin g et ap (Il997h and iBohringer et all (12002) . Right: Same as Fig. \T\ but for the REFLEX 
survey (iBohringer et al. 2004). One should note that the model predictions correspond to a 100% complete survey, while the actual 
survey is said to be at least 90% complete. 



the iJenkins et al.l ( 1200 ll) mass function), corresponding to 456 
clusters when considering 90% completeness and fully consis- 
tent with the observed number. Alternatively, considering our 
prediction as perfect would mean that the completeness of the 
survey is 88.2%. On the other hand, using erg, = 0.801 along 
with the Tinker et al. (2008) mass function leads to a total of 510 
clusters, hence 459 for a 90% completeness. A perfect predic- 
tion would, in this case, correspond to a completeness of 87.6%. 
Furthermore, as shown in Fig. |2] the overall shapes of the sim- 
ulated and observed redshift distributions are in both cases in 
very good agreement; note, however, that the completeness is 
not taken into account in this plot since the redshift of the miss- 
ing clusters is not known. 

Finally, concerning the SZ part of our model, we have 
already noted that the model provides an extremely good 
accounting of direct SZ obser vations (iMelin et al.l 1201 it 
iPlanck Collaboration et al.1 l2011dlla) . We hence consider the 
model to be a reasonable empirical representation of the current 
data and will now apply it to study the expected X-ray proper- 
ties of SZ-detected clusters. We emphasize that the model can 
be very easily adapted to new observational and theoretical con- 
texts: it is built around few analytical ingredients that can be 
trivially changed if needed, which is an important strength of 
this approach. 



3. SZ surveying and the Planck cluster catalog 

Like the X-ray emission, the SZ effect identifies clusters through 
the presence of the hot, ionized intracluster medium (ICM). The 
SZ effect, however, is less sensitive to gas substructure, being 
directly proportional to the thermal energy of the ICM, and re- 
mains unaffected by cosmic dimming. As a consequence, SZ 
cluster surveys efficiently find clusters at high redshift. 

The Planck satellite is the third generation space mission 
dedicated to the CMB, following CORE and WMAP. With its 
higher sensitivity, angular resolution and wide frequency cover- 
age, it is the first capable of finding large numbers of galaxy clus- 
ters through the SZ effect. Its all-sky catalog of clusters detected 
through the SZ effect (the Planck Cluster Catalog, or PCC) will 
be one of the primary and novel scientific products of the mis- 



sion. An early version of the catal og was recently published in 
iPlanck Collaboration et al.l ( 1201 l b, ESZ) based on high signal- 
to-noise detections in the first tens months of data. The Planck 
catalogs (ESZ and Legacy catalogs) are, remarkably, the first all- 
sky cluster catalog since those resulting from the ROSAT All- 
Sky Survey (RASS) of the early 1990s (Truemper 1992). They 
are expected to become a workhorse for many cluster and cos- 
mological studies. 

3.1. Properties of tlie Expected Pianck Cluster Catalog 
(EPCC) 

We use the matched-multifilter ("Herranz et al.l |20()2 |) detection 
algorithm described by Melin et al. (200g) to constitute the 
EPCC. Taking a set of sky maps at different frequencies, the 
algorithm simultaneously filters in both frequency and angular 
space to optimially extract objects with the thermal SZ spectral 
signature and the expected SZ angular profile. The filter is ap- 
plied with a range of angular scales to produce a set of candi- 
dates at each scale, which are subsequently merged into a single 
catalog. Construction of the filter requires adoption of a spatial 
template, for which we use the modified NEW pressure profile 
described in Eq. ([3]). In essence, we assume that the filter tem- 
plate exactly matches the true cluster profile, which should be 
borne in mind as an idealized situation. 

We apply our filter to simulations of the Planck data set 
based on an early version of the Planck Sky Model (PSIVQ)- 
These simulations provide a set of frequency maps including 
primary CMB temperature anisotropics in a WMAP-5 only cos- 
mology; Galactic synchrotron, free-free, thermal and spinning 
dust emission; extragalactic point so urces; sca n-modulated in- 
strument noise and beams (The Plan ck Consor tia 2005). Note, 
however, that we do not include any clusters in the simulations: 
by running our filter over these simulated maps, we seek only to 
quantify the total filter noise due to all instrumental and astro- 
physical sources. 

This gives us our noise threshold at each filter scale as a func- 
tion of position on the sky. With this information and our SZ 
scaling law, we then calculate the detectable cluster mass at a 



http://www.apc.univ-paris7.fr/ delabrou/PSM/psm.html 
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100001 




Fig. 3. Predicted cumulative redshift distribution of the EPCC 
for a nominal mission length (14 months) and threshold S /N > 
3. The beige band illustrates modeling uncertainty due to crfi and 
the choi ce of the m ass function. The solid lines are derived using 
the Jenkins et alJ j2001) mass function while the dashed lines 
correspond to the lTinkeret alJ (l2008l) mass function. 



given signal-to-noise threshold as a functio n of redshift, over the 
sky and for each filter scale. By using the iJenkins et al.l (1200 lb 
and Tin ker et aP P008) mass functions in our adopted cosmol- 
ogy, we find the EPCC distribution in mass and redshift. Figure[3] 
shows the resulting cumulative redshift distribution for S /N > 3 
at Galactic latitudes \b\ > 15° for a nominal 14-month mission, 
comprising two full-sky surveys, for three erg values that illus- 
trate the uncertainties on the determination of its value. It should 
be noted though, as shown above, that the cases that correspond 
best to the local counts are erg - 0.75 for th e lJenkins et al.l (1200 lb 
and (Tg = 0.80 for the [Tinker et al.l (l2008l) . Our results here are 
consistent with the more extensive study of a dozen detection al- 
gorithms run on simulated Planck data and presented in Melin 
et al. (2012, in preparation). 

3.2. Comparison with fhe RASS 

Since the Planck survey will be the first all-sky cluster survey 
since the RASS, it is of great interest to compare the EPCC 
to the "ROSAT Cluster Catalog", by which w e mean the sum 
of the REFLEX and NORAS catalogs (BohringereLaDHooS, 
l2004l) . This is done in Fig. HI which shows the measured ROSAT 
Cluster Catalog redshift distribution and the EPCC distribution 
for our three values of erg . From this comparison we see, as could 
be expected, that many clusters are common to both surveys. 
One of the first products of the Planck Cluster Survey will there- 
fore be the detailed study of the relation between SZ and X-ray 
properties afforded by the joint catalog. Such a study will fur- 
thermore benefit from the fact that the RASS clusters are gen- 
erally well-studied, allowing a number of important scaling re- 
lations to be derived. This will lead, for example, to additional 
tests of our empirical model and eventual modifications, all pro- 
viding useful insight into cluster physics. 

Moreover, prior knowledge of the X-ray-derived character- 
istics of these joint clusters will enable more precise recovery of 
their SZ properties. This is notably the case for cluster angular 







Fig. 4. Comparison of the ROSAT cluster catalog 
(REFLEX+NORAS; black) and EPCC differential redshift 
distributions, shown for different values of erg. The solid lines 
are derived using the Jenki ns et al] (1200 ll) mass f unction while 
the dashed lines correspond to the iTinker et al.l (|2008) mass 
function. The EPCC is expected to significantly extend the 
depth of all-sky cluster surveys. 



size, which is difficult to recover from the Planck data alone due 
to the size of the effecive SZ beam. This difficulty translates into 
noticably larger SZ f lux (Ysz) uncerta inties than if cluster size 
were known a priori (" Melin et al]|200Q) (also, Melin et al. 2010, 
in preparation). Clusters in the joint Planck-ROSAil catalog will 
therefore have significantly more precise SZ flux measurements 
thanks to the more precise X-ray information on cluster extent. 

More surprisingly, perhaps, is the fact that some RASS clus- 
ters are not seen in the EPCC. These tend to be mostly local and 
extended cl usters that are "resolved -out" by the Planck survey. 
As noted in lMelin etaP (I2005ll2006t) . the Planck selection func- 
tion depends on both cluster flux and extent. The fact that the 
Planck selection curve cuts through the observed RASS clusters 
is extremely important: it will be invaluable in constraining the 
PCC selection function by providing well-studied clusters that 
fall outside of the PCC and thereby improve our understanding 
of the selection criteria. This is a fortunate and unusual situation 
that we will fully exploit when constructing the PCC. 

Finally, as can be seen, we expect many new clusters not 
seen in the RASS. The PCC should greatly extend the redshift 
reach of all-sky cluster surveys over that of the RASS. The next 
section is dedicated to these new Planck clusters. 



3.3. XMM-Newton follow-up of the new Planck clusters 

We define new Planck clusters as those with a predicted X- 
ray flux in the ROSAT [0.1,2.4]-keV band less than /x = 3 x 
10"'^ ergs/s/cm^, which corresponds to the REFLEX detection 
limit and the deepest detection limit of the NORAS survey; we 
consider that all the other clusters are contained in the RASS cat- 
alogs. We then use the model to predict the X-ray properties of 
these new clusters in the XMM-Newton [0.5,2]-keV band. Fig. 
|5] shows the contours of isoflux within [0.15 - 1.0] x rsoo in the 
XMM-Newton band over the (z, T) plane for these new clusters. 
The key point is that we expect a substantial fraction of these 
new clusters to be hot (or equivalently, massive) and distant; 
in other words, those clusters that are the most useful cosmo- 
logical probes. Table |2] lists quantitative details on the distribu- 
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Table 2. Comparison between the predicted numbers of clusters in different subsets in the EPCC for different 0-^1, values. 



Characteristics 



Predicted number of clusters 















Jenkins 






Tinker 






(T, = 0.75 


cTg = 0.80 


0-8 = 0.85 


o-g = 0.75 


0-8 = 0.80 


<T8 = 0.85 


All 










1229 


2140 
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1316 


2 277 
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60 
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Fig. 5. X-ray characteristics of new EPCC clusters within [0.15 - 1.0] xrsoo- Left: Contours of iso-flux in the XMM-Afewfon [0.5-2]- 
keV band plotted in the (z, T) plane. Right: X-ray flux distribution of new high-redshift (z e [0. 8 - 1]) clusters for erg = 0.75 (blue), 
(Tg - 0.80 (green) and erg - 0.85 (red). The sohd lines correspo nd to results obtained using the lJenkins et al.l (12001 1) mass function 
while the dashed line was obtained using the [Tinker et al.l (l2008h mass function. 



tion of these clusters. It particularly focusses on hot (i.e. with 
T > 6 keV) and distant (i.e. with z > 0.5) clusters as they fall 
in the most poorly filled region of the (T, z) plane for known 
clusters. Moreover, these clusters are predicted to be very lu- 
minous, with X-ray flux in the XMM-Newton band greater than 
lO"'-' ergs/s/cm^ out to redshifts of order unity. 

The following expression provides a useful estimate of the 
exposure time texp (in ks) needed to measure a global tempera- 
ture to ~ 10% given its X-ray flux in the XMM bancfl: 



texp — JJ X 



/. 



10"'^ ergs /s /cm 



(10) 



We see that an exposure time of 55 ks is sufficient to obtain a 
temperature measurement for any cluster in the EPCC. Fig. |5] 
shows, as illustration, the histrogram of the number of new clus- 
ters with 0.8 > z> 1 .0 in the EPCC as a function of flux in the 
XMM-Newton band within [0.15-1 .0] x rsoo. All of the clusters 
have fluxes lying in the decade below the RASS detection limit 
and are thus easily observable with XMM-Newton or Chandra. 

In order to get an idea of the capabilities of XMM-Newton 
to observe such clusters based on actual measurements, we refer 
to the observations made of two of the four known clusters with 



M. Arnaud, private communication 



r > 6 keV and 0.8 > z > 1.0 and consider the results obtained 
and required exposure times: 



MS 1054-0321 JGioiaetal .1120041) : The XMM-A^ewfon obser- 
vations yield temperature, redshift (through observation of the 
iron line in the X-ray spectrum) and flux estimates of T - 



7.2t^ keV, 



0.847 



+0.057 
-0.040 



and /) 



X [0.5,2] keV 



(1.9 + 0.09) X 



10"'^ ergs/s/cm^, respectively; and they also provide a detailed 
image showing several components to the cluster's structure. 
The effective exposure time needed was ~25 ks. 



RX J1 226.9+3332 (Maughan et al. 2007): The same quanti- 
ties were measured in this case: T - 10. 4 + 0.6 keV, z - 0.89 an d 
/x [0.5,2] kev = 3.27 X 10-1^^ ergs/s/cm^ (IVikhlinin et al.ll2009al) . 
Moreover, the observation of the temperature profile, combined 
with the assumption of hydrostatic equilibrium, enabled a 
determination of the total mass and gas mass profiles, and thus 
an estimate of the cluster mass: M500 = 5.2;^Qg x IQ^*Mq. In this 
case, the exposure time was ~70 ks. 



These clusters are good examples of the kind of objects we 
find in the EPCC. We conclude therefore that 25-50 ks expo- 
sures is sufficient to obtain temperature measurements to ~ 10% 
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for any cluster in the EPCC, even out to z ~ 1 . It is even suf- 
ficient, in some cases, to obtain mass estimates. Detailed X-ray 
studies of a large fraction of the new Planck clusters is hence 
feasible and would significantly advance our understanding of 
cluster structure at intermediate to high redshifts. 

4. Conclusions 

The Planck SiZ survey will be the first all-sky cluster s urvey since 
the w orkhorse RASS dating from the early 1990s (iTruempen 
Il992h . The Planck Cluster Catalog (PCC) will thus satisfy the 
long-standing need for a deeper all-sky cluster survey. We have 
constructed an empirical model for cluster SZ and X-ray prop- 
erties that incorporates the latest results from detailed X-ray and 
SZ studies. Based on this model, we expect the PCC to extend 
the RASS to redshifts of order unity, as shown in Figs. [3] and |4] 
and to contain a high fraction of hot, massive and luminous clus- 
ters. These objects are the most useful for cosmological studies, 
because their properties (e.g., abundance) are the most sensitive 
to the expansion rate and their energetics is dominated by grav- 
ity, and less affected by feedback than are less massive objects. 

The expected X-ray properties of the PCC are given in 
Table |2] We see a significant number of new clusters, i.e., those 
not already contained in the RASS cluster catalogs (REFLEX 
and NORAS). Most should have temperatures T > 5 - 6 keV 
and have luminosities falling in the decade just below the RASS 
flux limit. Interestingly, they are easy to see with current X-ray 
instruments, if one knows were to look. And this is our key point: 
these useful clusters are rare and hence can only be found by an 
all-sky survey, such as the Planck survey. 

Our primary conclusion is that extensive X-ray follow-up of 
the new Planck clusters is feasible. Based on our cluster model, 
we estimate, for example, that XMM-Newton could measure the 
temperature of any newly-discovered Planck cluster to 1 0% with 
only a modest 50 ks exposure. This opens the path to impor- 
tant science with large follow-up programs on XMM-Newton and 
Chandra. A ~ 5 Ms program on XMM-Newton, for instance, 
could measure profiles and temperatures for ~ 100 Planck clus- 
ters at z > 0.6, and masses for a large subset. This would be 
a significant increase in sample size and lead to important ad- 
vances in the use of clusters as cosmological probes, as well 
as in our understanding of cluster physics at higher redshifts. 
While ambitious, this type of program is feasible and falls nat- 
urally under the category of Very Large Programmes envisaged 
for XMM-Newton after 2010. 

In conclusion, we expect the PCC to be unique for its ability 
to find rare, massive clusters out to high redshifts, and hence 
become a workhorse cluster catalog for many types of de- 
tailed cluster studies. In particular, we expect the PCC to ini- 
tiate important X-ray follow-up programs on XMM-Newton and 
Chandra. 

Acknowledgements. We are very appreciative of helpful conversations with M. 
Amaud and G. Pratt concerning cluster scaHng relations. We are also thankful 
for the very useful comments made by the referee that helped improve the quaUty 
of this paper A. Chamballu and J.G. Bartlett also thank the The Johns Hopkins 
University Department of Physics and Astronomy for their hospitahty during 
part of this work. A portion of the research described in this paper was carried 
out at the Jet Propulsion Laboratory, California Institute of Technology, under a 
contract with the National Aeronautics and Space Administration. 



Basu, K., Zhang, Y.-Y., Sommer, M. W., et al. 2010, A&A, 519, A29+ 
Birkinshaw, M. 1999, Phys. Rep., 310, 97 

Bohringer, H., Colhns, C. A., Guzzo, L., et al. 2002, ApJ, 566, 93 
Bohringer, H., Schuecker, P, Guzzo, L., et al. 2004, A&A, 425, 367 
Bohringer, H., Schuecker, P, Pratt, G. W., et al. 2007, A&A, 469, 363 
Bohringer, H., Voges, W., Huchra, J. P, et al. 2000, ApJS, 129, 435 
Burenin, R. A., Vikhlinin, A., Homstrup, A., et al. 2007, ApJS, 172, 561 
Burke, D. J., Collins, C. A., Sharpies, R. M., Romer, A. K., & Nichol, R. C. 

2003, MNRAS, 341, 1093 
Cole, S., Percival, W. J., Peacock, J. A., et al. 2005, MNRAS, 362, 505 
Dalai, N., White, M., Bond, J. R., & Shirokov, A. 2008, ArXiv e-prints, 803 
Davis, M., Efstathiou, G., Frenk, C. S., & White, S. D. M. 1985, ApJ, 292, 371 
Ebeling, H., Edge, A. C, Fabian, A. C, et al. 1997, ApJ, 479, L101 + 
Ebeling, H., Edge, A. C, & Henry, J. P 2001, ApJ, 553, 668 
Galh, S., Barflett, J. G., & Melchiorri, A. 2012, ArXiv e-prints 
Gao, L. & White, S. D. M. 2007, MNRAS, 377, L5 
Gioia, I. M., Braito, V., Branchesi, M., et al. 2004, A&A, 419, 517 
Halverson, N. W., Lanting, T., Ade, R A. R., et al. 2009, ApJ, 701, 42 
Herranz, D., Sanz, J. L., Hobson, M. P, et al. 2002, MNRAS, 336, 1057 
Hincks, A. D., Acquaviva, V., Ade, P A. R., et al. 2010, ApJS, 191, 423 
Homer, D. J., Perlman, E. S., Ebeling, H., et al. 2008, ApJS, 176, 374 
Jenkins, A., Frenk, C. S., White, S. D. M., et al. 2001, MNRAS, 321, 372 
Juszkiewicz, R., Feldman, H. A., Fry, J. N., & Jaffe, A. H. 2010, J. Cosmology 

Astropart. Phys., 2, 21 
Kaiser, N. 1986, MNRAS, 222, 323 

Komatsu, E., Smith, K. M., Dunkley, J., et al. 2011, ApJS, 192, 18 
Kotov, O. & Vikhlinin, A. 2006, ApJ, 641, 752 
Larson, D., Dunkley, J., Hinshaw, G., et al. 201 1, ApJS, 192, 16 
Lueker, M., Reichardt, C. L., Schaifer, K. K., et al. 2010, ApJ, 719, 1045 
Maughan, B. J., Jones, C, Jones, L. R., & Van Speybroeck, L. 2007, ApJ, 659, 

1125 
Melin, J.-B., Bardett, J. G., & Delabrouille, J. 2005, A&A, 429, 417 
Melin, J.-B., Bardett, J. G., & Delabrouille, J. 2006, A&A, 459, 341 
Melin, J.-B., Bartlett, J. G., Delabrouille, J., et al. 2011, A&A, 525, A139+ 
Menanteau, F, Gonzalez, J., Juin, J., et al. 2010, ApJ, 723, 1523 
Mohr, J. J., Mathiesen, B., & Evrard, A. E. 1999, ApJ, 517, 627 
Muchovej, S., Mroczkowski, T., Carlstrom, J. E., et al. 2007, ApJ, 663, 708 
Nagai, D., Kravtsov, A. V, & Vikhlinin, A. 2007, ApJ, 668, 1 
Navarro, J. E, Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 
Neumann, D. M. & Amaud, M. 2001, A&A, 373, L33 
Pacaud, F, Pierre, M., Adami, C, et al. 2007, MNRAS, 382, 1289 
Perlman, E. S., Horner, D. J., Jones, L. R., et al. 2002, ApJS, 140, 265 
Pifl'aretti, R., Amaud, M., Pratt, G. W., Pointecouteau, E., & Melin, J. . 2010, 

ArXiv e-prints 
Piifaretd, R. & Valdarnini, R. 2008, A&A, 491, 71 
Plagge, T, Benson, B. A., Ade, P A. R., et al. 2010, ApJ, 716, 1 1 18 
Planck Collaboration, Ade, P A. R., Aghanim, N., et al. 2011a, A&A, 536, Al 
Planck Collaboration, Ade, P A. R., Aghanim, N., et al. 2011b, A&A, 536, A8 
Planck Collaboration, Ade, P A. R., Aghanim, N., et al. 201 Ic, A&A, 536, Al 1 
Planck Collaboration, Aghanim, N., Amaud, M., et al. 201 Id, A&A, 536, AlO 
Pratt, G. W., Bohringer, H., Croston, J. H., et al. 2007, A&A, 461, 71 
Pratt, G. W., Croston, J. H., Amaud, M., & Bohringer, H. 2009, A&A, 498, 361 
Reichardt, C. L., Ade, P A. R., Bock, J. J., et al. 2009, ApJ, 694, 1200 
Romer, A. K., Nichol, R. C, Holden, B. P, et al. 2000, ApJS, 126, 209 
Staniszewski, Z., Ade, P A. R., Aird, K. A., et al. 2009, ApJ, 701, 32 
Sunyaev, R. A. & Zeldovich, Y B. 1970, Comments on Astrophysics and Space 

Physics, 2, 66 
Sunyaev, R. A. & Zeldovich, Y B. 1972, Comments on Astrophysics and Space 

Physics, 4, 173 
The Planck Consortia. 2005, The Scientific Programme of Planck ("Blue Book") 

(ESA Publication Division) 
Tinker, J., Kravtsov, A. V, Klypin, A., et al. 2008, ApJ, 688, 709 
Traemper, J. 1992, QJRAS, 33, 165 

Tytler, D., Kirkman, D., O'Meara, J. M., et al. 2004, ApJ, 617, 1 
Udomprasert, P S., Mason, B. S., Readhead, A. C. S., & Pearson, T. J. 2004, 

ApJ, 615, 63 
Vanderlinde, K., Crawford, T. M., de Haan, T., et al. 2010, ApJ, 722, 1180 
Vikhlinin, A., Burenin, R. A., Ebehng, H., et al. 2009a, ApJ, 692, 1033 
Vikhlinin, A., Kravtsov, A., Forman, W., et al. 2006, ApJ, 640, 691 
Vikhlinin, A., Kravtsov, A. V., Burenin, R. A., et al. 2009b, ApJ, 692, 1060 
Wu, J., Ho, P T. P, Huang, C, et al. 2009, ApJ, 694, 1619 
Zwart, J. T. L., Barker, R. W., Biddulph, P, et al. 2008, MNRAS, 391, 1545 



References 



Amaud, M., 
Amaud, M.. 



Pointecouteau, E., & Pratt, G. W. 2005, A&A, 441, 893 
Pratt, G. W., Pifl'aretti, R., et al. 2010, A&A, 517, A92+ 



